## Not run:
# suppressMessages(library(cubfits, quietly = TRUE))
#
# seq.data <- read.seq(get.expath("seq_200.fasta"))
# phi.df <- read.phi.df(get.expath("phi_200.tsv"))
# aa.names <- c("A", "C", "D")
#
# # Read in from FASTA file.
# seq.string <- convert.seq.data.to.string(seq.data)
# reu13.df <- gen.reu13.df(seq.string, phi.df, aa.names)
# reu13.list.new <- gen.reu13.list(seq.string, aa.names)
# y <- gen.y(seq.string, aa.names)
# n <- gen.n(seq.string, aa.names)
# scuo <- gen.scuo(seq.string, aa.names)
#
# # Convert to list format.
# reu13.list <- convert.reu13.df.to.list(reu13.df)
# y.list <- convert.y.to.list(y)
# n.list <- convert.n.to.list(n)
# ## End(Not run)
Run the code above in your browser using DataLab